arXiv:1508.02257vl [q-bio.NC] 10 Aug 2015 


Pre-stimulus oscillatory brain states and cognition: a theoretical approach 


Rakesh Sengupta 

Department of Computer Science and Electrical Engineering, 

Lassonde School of Engineering, York University, Canada 
(Dated: August 11, 2015) 

Spontaneous oscillations measured by local field potentials, electroencephalograms and magne¬ 
toencephalograms exhibits variety of oscillations spanning frrequency band (lHz-lOOHz) in animals 
and humans. Both instantaneous power and phase of these ongoing oscillations have commonly been 
observed to correlate with peristimulus processing in animals and humans. However, despite of nu¬ 
merous attempts it is not clear whether the same mechanisms can give rise to a range of oscillations 
as observed in vivo during resting state spontaneous oscillatory activity of the brain. In this paper 
I invetigate the spontaneous activity in the cortex. The paper attempts to establish analytically 
the conjecture that under certain conditions, a neural assembly can give rise to outputs that can be 
characterized by generalized oscillatory functions. It is possible to validate the analytical predictions 
with a neural mass model to show what the characteristic frequencies for such a brain state should 
be if they have to respond to external stimuli though it is not explored directly in the paper. In 
this paper we have attempted to show how an oscillatory dynamics might arise from a combina¬ 
tion of a feed-forward and recurrent neural assembly. Following that we have shown how naturally 
some of the EEG and MEG band activities in the pre-stimulus a (~ 10 Hz) can be explained from 
the resulting neural dynamics operating on a limited capacity cognitive systems. This provides a 
very important clue regarding how pre-stimulus brain oscillatory dynamics generates a window to 
consciousness. 


AUTHOR SUMMARY 

There is an inherent advantage to think of brain 
states following oscillatory dynamics. Recent works mag¬ 
netoencephalography and electroencephalography even 
contend that certain frequency ranges like a can be useful 
in understanding pre-stimulus brain states that precede 
conscious perception of stimuli. In this paper we attempt 
to show the theoretical possibility of thinking of oscilla¬ 
tory dynamics as a property of neural dynamics rather 
than just arising from properties of Fourier transform of 
signals. We have characterized the criterion for neural 
assemblies showing oscillatory dynamics and then pro¬ 
ceeded to show that theoretically there is justification in 
looking at the frequency bands like 9 and a. 


I. INTRODUCTION 

It has always been a standard practice to break down 
the brain signals obtained in electroencephalography 
(EEG) or magnetoencephalography (MEG) in terms of 
their frequency components obeying the laws of Fourier 
decomposition. It is even commonplace to compartmen¬ 
talise the frequencies into neat bunches like 9 (4-8 Hz), 
a (8-12 Hz), etc. However, the important question is 
whether these bands are just a matter of convenience or 
is there an underlying reality to the oscillatory model for 
the brain. Some recent work has tried to tie the oscilla¬ 
tory framework to underlying spiking model for neurones. 

S3 ■ This question can be solved only by looking at 
mathematical properties of neural codes. 

In recent years the understanding of neural codes has 
provided us with insights that go beyond the concepts of 


rate coding and it is increasingly more commonplace to 
speak of temporal codes that use spike timing and phase 
information in order to transmit and process information 
reliably S0SM1. However, most of these attempts 
are to look at neural codes after the presentation of stim¬ 
ulus. Recently, some studies have looked at pre-stimulius 
brain states in MEG and EEG based studies and have 
found that it is possible to predict conscious detection 
of stimuli based on pre-stimulus oscillatory brain activ¬ 
ity slim. For instance in case of near threshold 
stimuli some researchers have found the pre-stimulus a 
frequency band modulation to be important pbl ]. These 
attempts have drawn a large amount of interest, but have 
revealed little towards a theoretical or physiological un¬ 
derstanding of such phenomena. In the current work we 
have tried to start from a minimal number of assump¬ 
tions regarding a neuronal assembly and tried to show 
how it is possible to theoretically derive such an oscilla¬ 
tory dynamics. Moreover, we tried to capture the inher¬ 
ent constraints of the neural dynamics that can lead to 
a frequency region of interest in the a band. 


II. METHODS 

Conjecture 1: If a neural assembly S of N neurons con¬ 
sists of both feed-forward and recurrent connections 
with a finite variable bound refractory period r be¬ 
tween the feed-forward and recurrent connections, 
the equilibrium solution for the assembly can be 
characterized by a class of functions of the form 

A e i(uit-8) . 

Justification:: If a neural assembly comprises of ?zifeed- 
forward neurons and ri 2 recurrent neurons {n\ + 
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n 2 < N), then the general activations of the pools 
of neurons will be given by 

xi = ^-} ni = sgn(£ T £xi - 9) (1) 

X 2 = -^-}n 2 = -X2 + F.x 2 + r?I Xl (2) 

where £ is the weight matrix for the feed-forward 
connections, 9 is the threshold, rj is a scaling mul¬ 
tiplier, F is the non-linear transfer function and 
I Xl is the input coming from the feed-forward net- 
wor ME Now, for the network to be useful there 
needs to be sustained activity (provided by the re¬ 
current network, see EH) as well as the ability to 
restart the network dynamics. The latter is pro¬ 
vided by assuming that a smoothing regularizer is 
provided which varies inversely with the network 
output EH (to smooth the network against noisy 
perturbations), as well as a global decay parame¬ 
ter. In such a network, the overall output function 
< xi + x 2 > (average expectation value of the net¬ 
work output constructed in the line of mean activa¬ 
tion in varies between 0 and < x 2 >+/, where 
/ is the average normalized input to the network. 
This alone is sufficient to show the possibility of os¬ 
cillatory solutions. However, to be more rigorous, 
let us look at the second order time evolution for 
the feed-forward and recurrent populations. From 
Eq. U we get. 

xr = ~£r}rn = ±2<5(£ T £xi - 9) (5) 

Here 6(x) is the Dirac delta function. From Eq. [2] 
we have, 

dx 2 = -dx 2 + d{ F.x 2 ) 

= — dx 2 + (dF).x 2 + F ,dx 2 


If F is a smooth transfer function then the third 
term in the right hand side of Eq. [6] becomes neg¬ 
ligible. Thus we have 

x 2 = —x 2 + F.x 2 (7) 

x 2 = (l-F)x 2 +Fx 2 + 0(I Xl ) (8) 

For the recurrent network, near equilibrium (here 
we assume absence of input because of the finite 
refractory period of r), the nonlinear operators F 
and F become quasi-linear operators, and thus we 
have 

X 2 = ^}n 2 = -||F-(1+F)||x 2 (9) 

This has a form of simple eigenvalue problem. The 
network will have oscillatory solution if the opera¬ 
tor || F — (1 + F) || is positive Hermitian. A suit¬ 
able example is the leaky accumulator network de¬ 
scribed in (lll | where the matrix representing the 
operator becomes real symmetric. Thus near equi¬ 
librium, considering F is a quasi-linear approxima¬ 
tion < x 2 >, the mean state of the neural assembly 
near equilibrium / steady state will have a peri¬ 
odic solution £„ = Ae^ Unt ~ e ^ with xi contribut¬ 
ing to the phase 9. The full neural field dynamics 
if>(x, t ) can be constructed along the lines of Q with 
both spatial and temporal components taken into 
account as the following 

= E £ n (t) exp(inkx) (10) 

Thus the predictions from this analytical franework 
can be easily ported to neural mass models as well. 

Corollary: An assembly of recurrent neurons with 
complex-valued inputs and outputs, is formally 
equivalent to a neural assembly of independent 
feed-forward and recurrent neurons. 


1 If we consider a general recurrent shunting networks with dy¬ 
namics given by 

Xi — AiXi “t“ (Bi Xi)(Ii —h B(Xi'j'( 

(Xi + Ci) I Ji ) ' wjjS(Xj') 

V 

where I and J are excitatory and inhibitory inputs and S is 
a sigmoid function, we can transform it with symple variable 
change (yi = Xi + Ci) to the form 

b i {Vi) ~ E w oi S (.Vj ~ C 3 ) J ( 4 ) 




Justification: Let us start with Cohen-Grossberg gen¬ 
eralized networks of additive variant with a nonlin¬ 
ear activation function, like the network described 
by 111. If inputs and outputs to the net¬ 
work are given by a vector of complex numbers 
z{ixm} = x {lxm} + iy{ixm} for a network of m 
nodes, the network dynamics is governed by 

dz ■ 

= ~Zj + ciF(zj) - c 2 F(z k ) + Ij + noise (11) 

where F(z) = z/{\ + z). Ignoring noise for the time 
being, if we decompose the real and imaginary parts 
of EH we have 


where bi(yi) = 2- {AiCi - (A; + Ji)yi + (Bi + Ci - yi) [h + S(yi - Ci)]}. 
It can be shown that in absence of input and near equilibrium, 
second order dynamics of this network will be very similar to 
Eq. \S\ with an extra term of the order of y ?, which does not 
change the general conclusions of the paper. 


dz-j Zj + lzA 

= + Cl l + N 2 + 2R e(^) 


-*E 


Zk + \Zk \ 2 


k^j 1 + |~fc| + 2Re(^fc) 


+ I-i 


( 12 ) 
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And thus the separated real and imaginary parts 
yield two equations given by, 


dRe(zj) 


dt 

-*£ 


= -Re(zj 
Re(zfc) 


Ci 


Re(zj) 




1 + |~fe| + 2Re(2:fc) 


1 + kj'l +2Re(zj) 
-Re(/,) + 0(|z| 2 ) 


(13) 


= Y c jk lla ( z k) + (14) 

where c' fe are constants. Making suitable substitu¬ 
tions (xj f— 2(Re(^) + |zj| 2 /2), and yj 4— Im(^) 
), we have 

= ~Xj + c"F(xj) - 4' Y F ( x k) + !j 

dt fe (15) 

-(-constant terms 


■Jr - £ c 'jkyk +( 16 ) 

From Eq. [T5] and [TH] it is clear that a combination 
of recurrent and feedforward networks can function 
in a way to handle complex inputs and outputs and 
thus at least the reverse of Conjecture 1 is true in 
certain casefl 

Conjecture 2: If a neural assembly § of N neurons con¬ 
sists of both feed-forward and recurrent connections 
with a finite variable bound refractory period r be¬ 
tween the feed-forward and recurrent connections, 
the pre-stimulus brain states are determined by the 
delay between the regions represented by the neural 
assemblies. 

Justification: The conjecture 1 has some interesting 
consequences. Firstly, it allows for brain states 
to be defined in terms of a frequency or a distri¬ 
bution of frequencies. In a case of a state with a 
distribution of frequencies we can think of a charac¬ 
teristic frequency range representing the state, the 
characteristic frequency being the one with maxi¬ 
mal power. The conjecture also allows us to think 
of the brain states as phenomenal superpositions 


2 Thie above analysis also holds for a general sigmoid activation 

function. If <r{x) = | | )_ r , : , then from simple function approxi¬ 
mation of a general sigmoidal function f(x) = c /. a (x — a,) 

can be simply transformed in the complex domain as f(x) —> 

z Z +a i where cti = e“* (see Appendix C of f(J). Thus the 
results of conjecture 2 are quite general. 

3 Here we are using the term brain state to mean the state of the 
neural assembly involved in a particular task or function, not the 
entire brain. 


of oscillatory dynamics, allowing us to deal with 
problems such as stimulus related perturbations to 
brain states more efficiently. 

The prevailing additive idea of brain states in 
the neuroimaging literature needs no introduction. 
However, we will formally spell out the bare essen¬ 
tials. Before a test condition appears to a subject 
the family of brain states or the neural assembly 
§ in question ( S ) can be thought of as its resting 
state. In the test condition, a new perturbation 
comes from our experimental control (St)- Or we 
can write, S <r- S + St- 

Now if the resting state is imagined as a stand¬ 
ing wave, then we have from Eq. [10] S = 
Acos(kx)exp(iojot). Since resting state can be 
taken to be not very location specific we can write 
the rate of change of the state of neural assem¬ 
bly given by S, when a perturbative brain state St 
(generally due to oncoming stimulus) interacts with 
the current state to be given by 

— = *aexp(u;ot) + (17) 

From the results of the previous section, if assume 
the brain states and their perturbations to have 
solutions of the form we can write, 

ai iwe i(ajt " e) = iae^ + a 2 ico'e i( - Jt ~ 8 ^ (18) 

Here w is the pre-stimulus brain state frequency 
that tries to interact with the incoming perturba¬ 
tion characterized by u/. Separating the real and 
the imaginary parts and applying the constraint 
that the imaginary part must go to zero on the left 
and right hand side of Eq. [lH] 

aiw cos (ut — 9) = a cos(ujot) + a 2 w' cos (u't — 9±) (19) 

Now considering the case to ~ a/, we have from Eq. 
l22l (if S and St are in similar phase) we have for 
t -A 0, 


(ai — a 2 )uj cos(ujt — 9) = a 


Thus we have our constraint, 


- 1 < 


a 

(a i — a 2 )co 


< 1 


( 20 ) 


( 21 ) 


Thus Eq. HU shows that the resonant frequency of 
the pre-stimulus brain states in a region varies with 
a 0 . Interestingly, a is a spatial term dependant on 
the spatial connectivity as seen above. Thus 


4 Now if we take the imaginary part of the Eq. 1221 then we have 
(for the case u) ~ oj' if S and St are in similar phase) 
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ut - 9 ~ (25) 

O(ojt) ~ n7r (26) 

If S' and St are anti-phase then also we have 
tan(wf — Q) = 2± <oc and thus 0{u)t) ~ 7r.On 
the other hand for the condition w u/, we have 
cos (cot — 9 ) =0, and thus Eq. [24] follows. The 
condition a/ w is not very informative as it as¬ 
sumed very high frequency stimulus for low fre¬ 
quency brain states. Thus overall we have that 
a pre-stimulus oscillatory brain state characterized 
by the frequency w, will be informative about the 
stimulus if the Eq. [24] holds. Interestingly, most 
brain dynamics is a time limited phenomena, and 
the frequency w thus depends upon the time con¬ 
stant of the feedback connection (r). However these 
are related to the biophysical process rates. Thus 
this formulation gives a very natural way of looking 
at which frequencies might appear in the dominant. 
For instance, if r = 100 ms, then v = ^ = - would 
be of the order of 10 Hz. This idea is derived here in 
a non-trivial manner, i.e., just having a refractory 
period of 100 ms does not guarantee oscillatory so¬ 
lutions. Here we have shown how it is possible to 
have oscillations in the a band can be thought to 
arise from internal brain dynamics. 

Thus there are two consequences of the analytical 
framework. It shows how a general oscillatory activity 


can be generated in the particular brain region having 
both feedforward and recurrent connections. Secondly 
it connects the general oscillatory signals in the brain 
that arises from connected regions to be dependant upon 
the delay in the network connectivity arising from spatial 
factors. 

III. RESULTS 

In the above, we have shown the possibility of oscil¬ 
latory brain states under certain equilibrium conditions 
for a neural assemble consisting of both feed-forward and 
recurrent connections. We also showed how the oscilla¬ 
tory characterization of the states in the brain or neural 
assembly leads quite naturally to understand the evolu¬ 
tion of brain states in terms of superposition of states. 
This formalism ultimately leads to the result that brain 
states in the pre-stimulus domain can interact or resonate 
with the post-stimulus perturbations if the frequencies 
of the states are of the order of value or its multiples. 
However, it still does not say how far in the past the 
pre-stimulus state needs to be to have the desired effect. 
These have to be further explored by implementing neu¬ 
ronal mass models to verify the predictions. 
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a. 2 ) sin(t ot — 6) = 0 

(22) 

cot — 0 ^ Tin 

(23) 

O(cot) ~ nn 

(24) 


If S and St are anti-phase then also we have tan (uit—0) = <oc 
and thus 0(ut) ~ 7r.On the other hand for the condition uj '^> N, 
we have sin(ojt — 6) = 0, and thus Eq. 1241 follows. The condi¬ 
tion uS uj is not very informative as it assumed very high 
frequency stimulus for low frequency brain states. Thus overall 
we have that a pre-stimulus oscillatory brain state characterized 
by the frequency lj, will be informative about the stimulus if 
the Eq. 12II holds. Interestingly, most brain dynamics is a time 
limited phenomena, and the frequency lu thus depends upon the 
delay between the connections as shown in 1211 However these 
delays are related to the biophysical process rates. Thus this for¬ 
mulation gives a very natural way of looking at which frequencies 
might appear in the dominant. For instance, if the delay r = 100 
ms, then v = -£- = — would be of the order of 10 Hz. This idea 
is derived here in a non-trivial manner, i.e., just having a refrac¬ 
tory period of 100 ms does not guarantee oscillatory solutions. 
Here we have shown how it is possible to have oscillations in the 
a band can be thought to arise from internal brain dynamics. 
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